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Abstract 

Proportional-Integral-Derivative (PID) controllers are used in many of the Industries for various 
process control applications. PID controller yields a long settling time and overshoot which is not good for 
the process control applications. PID is not suitable for many of the complex process control applications. 
This research paper is about developing a better type of controller, known as MPC (Model Predictive 
Control). The aim of the paper is to design MPC and PID for a pasteurization process. In this manuscript 
comparison of PID controller with MPC is made and the responses are presented. MPC is an advanced 
control strategy that uses the internal dynamic model of the process and a history of past control moves 
and a combination of many different technologies to predict the future plant output. The dynamics of the 
pasteurization process was estimated by using system identification from the experimental data. The 
quality of different model structures was checked using best fit with data validation, residual and stability 
analysis. Auto-regressive with exogenous input (ARX322) model was chosen as a model structure of the 
pasteurization process and fits about 80.37% with datavalidation. MPC and PID control strategies were 
designed using ARX322 model structure. The controller performance was compared based on settling 
time, percent of overshoot and stability analysis and the results are presented. 
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1. Introduction 

In a modern world the economic and quality issues become more and more important, 
efficient control systems have become indispensable. Therefore the process industries require 
more reliable, accurate, robust, efficient and flexible control systems for the operation of process 
plant. In order to fulfill the above requirements there is a continuing need for research on 
improved forms of control [1]. 

Control of temperature plays an important role in pasteurization plants. High 
temperature short time (HTST) is keeping milk or other food stuffs at 72 OC for 15 seconds in 
insulated holding tube. The pasteurization process consists of three stages like regeneration, 
heating and cooling sections. The crucial stage is heating process using heat exchanger to 
ensure unpasteurized product achieve desired pasteurization temperature before pass through 
holding tube and cooling sections. Prior to pasteurize milk sample, the equipment must have 
adequate controller to control the outlet temperature in order to maintain at standard value [2]. 

The proportional integral (PI) and proportional integral derivative (PID) controllers are 
widely used in many industrial control systems because of its simple structure. This controllers 
are designed without process constraints only use mathematical expression based on error from 
a set point. In these circumstances, conventional controllers (PI and PID) are no longer to 
provide adequate and achievable control performance over the whole operating range. Thus 
designing a controller considering the process constraints and optimize the control performance 
is essential [3]. 

Model Predictive Control also known as receding horizon control, is an advanced 
strategy for optimizing the performance of multivariable control systems. MPC generates control 
actions by optimizing an objective function repeatedly over a finite moving prediction horizon, 
within system constraints, and based on a model of the dynamic system to be controlled [4], 


Received November 7; 2016, Revised January 7, 2017; Accepted January 21, 2017 



Bulletin of EEI 


ISSN: 2302-9285 


■ 25 


2. Process Description 

The plant PCT23, manufactured by Armfield (UK), is a laboratory version of a real 
industrial pasteurization process, as shown in Figure 1. It consists of a bench-mounted process 
unit to which is connected a dedicated control console. An interface card DT2811 is used for 
monitoring and controlling the process through a computer [5]. 



Figure 1. Process flow diagram of pasteurization plant 


Here temperature T1 is the controlled variable and milk flow rate (N1), hot water flow 
rate (N2) and power are manipulated variables. 


3. Experimental Setup 
3.1. Input-Output Data 

The input-output data was generated by introducing step input in manipulated variables 
(N1, N2 and power input) and record temperature response (T1). Two input-output data was 
generated for model estimation and validation as shown in Figure 2. 

3.2.System identification 

The input-output data was analyzed by the System Identification toolbox in MATLAB. 
The continuous and discrete model structures were tried to select the model structure that have 
best fit with validation. Then the selected model structure is tested for residual analysis and 
pole-zero analysis to check the model stability. The continuous time (CT) model, Auto¬ 
regressive with exogenous input (ARX) model structures, ARMAX (auto regressive with moving 
average and exogenous (or extra) input model, Output-Error (OE) model and state space model 
structures were tried get best model structure in terms of best fit with validation data and model 
stability for further controller design [6]. 
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Figure 2. Systematic procedure for model identification process 


Best fit is calculated as: 

Validation data (y) is estimated data and is mean of validation data. After selection of 
best fit model structures model quality analysis like residual and pole zero location should be 
checked to select a nice and simple model for further controller design. The prediction error or 
residual is the key quantity. It is defined as the stability of a system can be easily inferred by 
examining the pole locations of the transfer function [7], 

3.3. Controller Design 

Controllers are basically employed in a closed loop control system. Closed loop control 
system is one that automatically changes the output based on the difference between the 
feedback signals to the input signal. Controller is an element used to produce manipulated 
variable from error variable, for control action [8-9]. 

a. PID controller 

The error signal e(t) is the difference between the reference input r(t) and desired output 
y(t). This error is manipulated by the PID controller to produce a command signal for the 
system. PID controller is described by the following transfer function in the continuous s- 
domain. 

b. Model Predictive Control 

The model predictive uses quadratic minimization problem. Subject to input and output 
constraints of the system, where is the set point, Q1 is output weight and Q2 is input weight 
as shown in Figure 3. The size of this minimization function and weight matrixes are depend 
on prediction and control horizon [10-11.] 


4. Results and Discussions 
4.1. Input-output data 

Two different experiments are done for model estimation and validation purposes at 
sampling time of 10 seconds. 
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Figure 3: Input-output data profile 


4.2. Model Structure selection 

The continuous time model fits the validation data better than the others. ARX422, 
ARMAX3202 and state space model pss2 are also fits the validation data fairly. The continuous 
time model doesn’t mean a good model rather the selection of model needs further analysis, as 
shown in Figure 4. 



Figure 4. Fits of different model structures 


4.3. Model Quality Analysis 
4.3.1. Residual Analysis 

For different model structures the auto corelation of residuals for the output (whitness 
test) and cross correlation of residuals with the input (independence test) were analyzed. From 
the graphs the horizontal scale is the number of lags, which is the time difference (in samples) 
between the signals at which the correlation is estimated. The yellow region on the plot 
represents the confidence interval of the corresponding estimates. Any fluctuations within the 
confidence interval are considered to be insignificant. A good model should have residuals 
uncorrelated with past inputs (independence test) and past outputs (whitness test). For poor 
models either auto and cross corelation residuals or two of residuals is out of the yelllow or 
confidence regionion as shown in Figure 5. 
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Figure 5. Residual analysis for different model structures 


Only auto correlation analysis of continuous time model analysis shown as some 
residuals is out of the confidence region. This means this model fail the residual analysis and 
the other three models still needs other model quality analysis to choose which model is better 
to represent the pasteurization process. 

4.3.2. Pole-Zero location 

Zeros and the poles are equivalent ways of describing the coefficients of a linear 
difference equation. Poles are associated with the output side of the difference equation, and 
zeros are associated with the input side of the equation. Generally for stable process all of the 
poles of discrete time model should be found within unit circle otherwise the model is not stable 
as shown in Figure 6-8. 



Figure 6. pole-zero location of ARX422 model 
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Figure 7. Pole-Zero Location of ARMAX3202 model 



Figure 8. Pole-Zero Location of State Space Model 


The ARX422 model is stable and there is a chance for pole-zero cancelation. Therefore 
we can reduce the model order. The ARMAX3202 also stable but the state space model is 
marginally stable. This model is rejected here. 

4.4. Model Reduction 

The ARX422 reduced to ARX322, still the fit is greater than ARMAX3202. Further 
reduction below this order deteriorates the fit. The reduced model (Figure 9) also passes the 
model quality analysis. Therefore the ARX322 model represents the pasteurization process. 
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Measured and simulated model output 



Figure 9. Reduced Model Fit 


The following equation is the converted discrete ARX model to continuous dynamic 
model for PCT23 pasteurization plant. The process has three inputs that come parallel and one 
output. The process dynamics is the third order process. 

Input delays are 10, 80 and 30 respectively for the three inputs (milk flow rate, hot water 
flow rate and power input). The arx322 model has Loss function of 0.0334283 and final 
prediction error (FPE) of 0.034922. 

4.5. PID Controller 

The milk flow rate is set to constant, because our pasteurization is high temperature 
short time pasteurization that takes only 15 seconds. In order to meet this requirement we have 
to set our milk flow rate at 326 ml/min. The hot water and the power input are let for PID 
controller (Figure 10) to adjust them until the set point (72 0C) is reached. 




Figure 10. Simulink block diagram of PID controller 
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The two PID controllers tuned by trial and error method until the set point is tracked 
efficiently by using Simulink software. Table 1 is the setting of the two PID controllers to get best 
set point tracking. 


Table 1. PID controller setting 


Setting 

PID controller 1 

PID controller 2 

Kc 

0.06 

0.3 

Ti 

0.12 

0.0009 

Td 

1 

0.9 


Based on this setting the output temperature response is shown in Figure 11. 



Figure 11. Output Temperature (T1) Response Using PID 


4.6. MPC Controller 

The prediction horizon of the system should be large enough to cover the settling time 
as shown in Figure 12. The prediction horizon of 200 samples gives better settling than the 
others prediction horizon intervals as shown in Figure 13. 
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Figure 12. Simulink block diagram for MPC design 
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Figure 13. Prediction Horizons Tuning for MPC 


The control horizon for different system should be different. It depends on the output 
signal of the system. In the most cases the control horizon should be large enough to get the 
reasonable stabilize output signal of the system. The long control horizon is required to improve 
the performance. Control interval of 95 is taken because beyond that it makes overshoot and 
becomes unstable, as shown in Figure 14. 



Figure 14. Controller Horizon Tunning for MPC 


The output weight is 1 and input weight is zero. For the above values of prediction and 
control horizons, input and output weights and plant constraints, the optimized output 
temperature response using MPC is shown in Figure 15. 
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Plant Output T1.°C 
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Figure 15. Output Temperature Response Using MPC 


4.7. Stability Analysis of MPC 

The stability analysis of MPC can be checked using pole-zero location of the closed 
loop response after using MPC. All the poles and zeros are existed inside unit circle for the 
three inputs to one output relations as shown in Figure 16. 



Figure 16. Closed loop pole-zero locations 


4.8. Comparison MPC and PID Controller 

The output temperature should be maintained at 72 0C which is the optimum 
temperature for high temperature short time (HTST) pasteurization process, but PID controller 
have overshot above this temperature, while MPC targets the set point exactly. The rise time of 
PID controller is better than MPC, but this may be due to high adjustment of power input out of 
the constraint region as shown in Figure 17. 
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Figure 17. Output pasteurization temperature response 


MPC can better set point tracking ability than PID controllers because this controller 
working in optimal condition by considering process constraints, that why there is no overshoot 
during MPC control action. Both controllers have the same settling time. MPC produces no 
overshoot while PID have higher overshoot. 


5. Conclusion 

From the above results and discussion we identified the best controller for 
pasteurization process to be MPC. Maintaining the temperature at a constant value is a critical 
issue in many of the Industries. MPC fulfills these type of difficulties by bringing the process 
variable to the desired set point as early as possible. MPC controller are more suitable or 
complex process control applications. MPC controller has 0% overshoot and minimum settling 
time as compared with PID controller that we used in this process. Consolidating the 
performance of all these responses MPC is the best controller for our pasteurization process. 
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